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Abstract. We consider positive radial solutions to the fast diffusion equation ut = Au™ on 
the hyperbolic space for N > 2, m € (nis, 1), — By radial we mean solutions 

depending only on the geodesic distance r from a given point o £ El'^. We investigate their fine 
asymptotics near the extinction time T in terms of a separable solution of the form V{r,t) — 
(1 - f/r)i/(i-'"V'''"(r), where V is the unique positive energy solution, radial w.r.t. o, to 
— AV^ = c y^l™- for a suitable c > 0, a semilinear elliptic problem thoroughly studied in |29| . [7]. 
We show that u converges to V in relative error, in the sense that Ijit^C-, i)/V"'(-, t) — ^ 
as f — ^ T~ . In particular the solution is bounded above and below, near the extinction time T, 
by multiples of (f — t/r)i/(i-'")e~('^~i)''/™. Solutions are smooth, and bounds on derivatives 
are given as well. In particular, sharp convergence results as t — >■ T~ are shown for spatial 
derivatives, again in the form of convergence in relative error. 



1. Introduction and preliminaries 

We analyze the asymptotic behaviour of solutions to the following Fast Diffusion Equation 
(FDE) on the A^-dimensional hyperbolic space (throughout the whole of this work we shall 
always assume N >2): 

jut = A (tt™) on X (0, T) 
|n = lio > on X {0} ' 

where A is the Riemannian Laplacian, m £ {nis, 1) and is the critical exponent 

A^- 2 



The initial datum uq is assumed to be radial in the sense that it depends only on the geodesic 
distance r from a given point o (which we shall indicate as d{x, o)), such a point being considered 
as fixed. Solutions to the FDE corresponding to radial data are of course radial as well for any 



fixed time. In (1.1) the parameter T = T(uo) denotes, for an appropriate class of data, the 
extinction time of the solution u, namely the smallest positive time t at which u{x, t) = 
identically. In fact, the results of [8] show that, in a class of Cartan-Hadamard manifolds 
which includes H^, such a time does exist finite for initial data which belong to L''(EI^), where 
q > max (1, N{1 — m)/2). We refer to such paper also for the relevant existence and uniqueness 



results provided there for solutions to (1.1) (see however an alternative approximation procedure 
sketched in Section [2]). 

It is well-known that the fine asymptotics of solutions to the FDE posed in the whole Eu- 
clidean space is governed by (suitable rescalings of) Barenblatt, or pseudo-Barenblatt, 
solutions. A huge literature on the topic has been produced in the last decade, and we 
limit ourself to quote the recent book |34j and, without any claim of completeness, the pa- 
pers [l5l[l3l|18l|ll|26l|2ni[l2l[3Ol[I71|5llS0 and references quoted therein. Notice 
that extinction in finite time holds, in this context and for appropriate class of data, only for 
m < {N — 2)/N (provided one takes, for simplicity, > 3). 

The situation on negatively curved manifolds is different from the Euclidean one since van- 
ishing of solutions in finite time often occurs not only for m close to 0, but for all m < 1. In 
fact, this is instead somewhat similar to what happens in the case of the homogeneous Dirichlet 
problem on bounded Euclidean domains il., which was deeply investigated in \21\ l23l [2l (TUl I32j . 
where at various levels of detail it is shown that the asymptotics of suitable classes of solutions 
can be discussed in terms of separable solutions of the form (1 — t/T)^/^^^''^'^ S^^'"^ (x) , S being a 
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positive solution to the elliptic problem —AS* = cS^^^ in fi, S" = on dQ, for an appropriate 
value of the constant c > (in principle depending on the initial datum). 

No positive solution to the above elliptic problem exists in the whole Euclidean space M^, 
but the fact the the bottom of the spectrum of —A on H'^ is strictly positive points towards 
existence of such a solution in . In fact, this result has been proved in [29] , More precisely, 
it is shown there that, given any c > and m £ {rus, 1), the equation 

-AV = cV^ onM^ (1.2) 

admits strictly positive solutions V belonging to the energy space H^m^) (what we call energy 
solutions) , which are necessarily radial with respect to some point o, the latter being therefore 
the only free parameter characterizing such solutions. Notice that solutions to (1.2) associated 



to different values of c are related by scaling, namely they are all multiples of the solution 
corresponding to c = 1. The asymptotics of y as r — ?■ oo has been studied as well in |29) and 
slightly improved in [7]: the main result is the existence of constants I = l{c,m,N) > such 
that, for all k G NU{0}, 



lim e 

r— ^-oo 



Actually, in [29] and |[7j, ( 1.3 ) is shown in to be valid for k 
by V allows to prove it for any k gN. 



(r) = (-l)'=(Ar_i)fc/. (1.3) 
0, 1 only, but the equation satisfied 



Infinitely many other positive solutions to (1.2) exist, but none of them belongs to H^{M^) 
and their behaviour as r — t- oo is polynomial (see again [7J). 

The asymptotics of a given solution u to ( |1.1[ ) starting from a radial initial datum uq is 
related to the energy solution V of (1.2) having the same pole a as uq and corresponding to a 



choice of c that depends on uq itself via the extinction time of u, namely the one that satisfies 

-AV=-, ^——V^ onM^. (1.4) 

(l-m)T ^ ' 

In fact, our main result is the following Theorem |1.1[ Formula (1.5) in it will be proved in 
Section [5] while formulae (1.6) and (1.7) will be proved in Section |6] 

Theorem 1.1 (convergence in relative error and convergence of derivatives). Let u be the 
solution to the fast diffusion equation (1.1) corresponding to a non-identically zero initial datum 
Uq > 0, which is radial w.r.t. o S M.^ and belongs to L'^{Il^) for some q > A^(l — m)/2 with 
q > 1. If T > is the extinction time of u and V is the unique positive energy solution, with 



pole 0, to the stationary elliptic problem (1.4), then 



lim 

Moreover, for all k £N there holds 



lim 

where it is understood that 



u{t) 



[l - y-, 



u{t) 



1 - 1,) 1-- 



0, 



(1.5) 



(1.6) 



C"=(R) 



for all regular functions (p. As a consequence, for all k £ N there exists a smooth function 
Fj.{r), having the property 

lim Fk{r) = 1 , 
1 — ^00 

such that 



lim 



0. 



(1.7) 
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Remark 1.2. Formulas (1.7), (1.3) (and identity (6.15) below) imply that, given k GN, for all 
e > there exist t^: G (0, T) and > such that 



1 - e < 



(1 



-(r) 



< 1 + e Vr >r, , Vt G [t„T) . 



(1.8) 



Notice that (1.7) bears some similarity with some of the results given, for the Euclidean case 



and in the range of m for which there is no extinction, in |28] . See also |27) for similar results 
for the Euclidean p-Laplacian driven evolution. 

The method proof of Theorem 1 1 . 1 1 and the known behaviour at infinity of V and its derivatives 
allow to show that the next Theorem |1.3| holds. In it, we shall first state a global Harnack 
principle, in the spirit of |21] and |11| . Secondly, we shall give upper and lower bounds on 
derivatives of the solution. In fact, (1.9) will follow from Propositions 3.1 and 4.1 (l.lO)-(l.ll) 
follow from the results given in Section [6] whereas (1.12) is an immediate consequence of (1.8), 
(LSl) and (felsl). 



Theorem 1.3 (global Harnack principle and bounds for derivatives). Let the assumptions of 
Theorem be valid. Then for all e > there exist positive constants ci = ci{uQ,m,N,e), 
C2 = C2{uo,m, N, e) such that the bound 



ci 1 



t \ 1-™ _N-1 ft 
7f] e ra < u{r, t) < C2 ( 1 - - 

holds true for all r >0, t G [e, T). Moreover, for all k gN there hold 

Qk 



1 

1 — m 



Qk 



Qj.k 



u{r, t) 



u{r, t) 



1 



1 

1 — m 



< G 



2,k 



'I 



T' 



1 

1 — m 



(1.9) 

[1.10) 
[1.11) 



Vr > , Vt G [e, T) 



for suitable positive constants Ci^k = Ci^k{uQ-,m,N,e), C2,k = C2.kiuo,m,N,£). In addition, 
for any k gN there exist i G (0,T), f > and a suitable positive constant C'^ = C'i^{uo,m, N) 
such that 



Qk 

Qj.k 



u{r, t) 



> Ci 



1 



1 

1 — m 



Vr>f, VtE [t,T). 



[1.12) 



Notice that our global Harnack principle (as well as the corresponding bounds for derivatives) 
is in the spirit of the one proved in the fundamental paper |i21j by DiBenedetto, Kwong and 
Vespri for bounded domains of M^, and of the corresponding results proved by Bonforte and 
Vazquez |11| for the fast diffusion equation in M^: in this latter case solutions can in fact be 
bounded above and below for all times by suitable Barenblatt solutions. Convergence in relative 
error was first discussed in ^33j for solutions to the fast diffusion equation in M^, the attractor 
in that case being still a Barenblatt solution. Later, Bonforte, Grillo and Vazquez showed in 
[10] that convergence in relative error to a separable solution occurs in the case of bounded 
domains, thus improving the results of |21) . 

It is worth pointing out that the techniques of proof of Theorems 1 1 . 1 1 and 1.3 can be used to 
capture the spatial behaviour of solutions, for any fixed t > 0, to the FDE also in the subcritical 
range m G (0,ms]. Indeed the following result holds (see Remarks |3.8[ [4^ and 6.1 for a sketch 
of proof). 

Theorem 1.4 (spatial behaviour for subcritical m). Let the assumptions of Theorem be 
valid, and suppose that m lies in the subcritical range (0,ms]. Then for any fixed t G (0,r) 
there exist positive constants ci = ci{t,uo,m, N), C2 = C2{t,Uo,m, N) such that the bound 
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holds true for all r >0. Moreover, for all k gN the bounds 

Qk 



u{r, t) 



Qk 



1.14) 
1.15) 



Vr > 0, vt e (o,r) 

hold true for suitable positive constants Ci^k = Ci^k{t,uo,m, N), C2,k = C2,k{tTUi),m, N). 

Some words have to be said about the assumption of radiality that we require on the initial 
data, which is related to two technical issues. Firstly, we need some a priori decay properties 
for the solution in order to exploit suitable barrier arguments. Such decay properties hold 
automatically for radial functions in the energy space but need not be valid for general solutions. 
In second place, it is not obvious that the solution (suitably rescaled in time) corresponding 
to a nonradial datum selects a unique limiting spatial profile V along subsequences (recall the 
degree of freedom given by a). This was proved in |23] in the Euclidean case (on bounded 
domains) but it is not known up to now in the present context. 

Remark 1.5. Keeping the fundamental hypothesis of radiality, our results hold in somewhat 



to 



more general geometric frameworks, but we preferred to state them in the case of 
avoid bothering the reader with heavier notation and technicalities. In fact one could consider 
Riemannian models (see |25| [3] as general references and |1J for the analysis of Lame-Emden- 
Fowler equations in such context) whose metric is defined, in spherical coordinates about a pole 
o, by ds"^ = dr'^ + iP^{r)de^, G G where V G ^^([o^oo)), V(0) = <(0) = 0, ^'(0) = 1, 

^l^'{r) > for every r > and limr_5>oo V''('")/V'(^) ^ (0, cxd). Notice that sectional curvatures at 
a point P tend, as the geodesic distance d{o, P) tends to oo, to a strictly positive constant. In 
such a kind of manifold a radial energy solution having the properties of the present solution 
V has been shown to exist in jl]. 



1.1. Preliminaries. As for the initial datum uq = Uo{r) in (1.1), in principle besides its 
nonnegativity we should also assume that it is bounded and such that ti™ G H^g_^{M^), where 



V radial: 



v\s) (sinh sr-Ms + / [v'{s)Y (sinh sr'Ms < oo 



(1.16) 

Notice that H^^^{Il^) coincides with the space of radial functions about o which belong to 
H^{Il'^). By energy solutions to (1.1) one should mean those starting from data uq as in (1.16) 
(in the nonradial context, those starting from uq : S H^{M^)), but in fact the results of 
[H] show that the solution u corresponding to an initial datum which fulfils the integrability 
conditions of Theorem O automatically satisfies n"'(-,e) e H^^^{Il^) n L~(H^) for all e > 0. 
This is stated in j8j for > 3, but it holds true when A = 2 as well because the methods 
of proof exploited in |8| rely only on the validity of a suitable Sobolev inequality in H^{M'^), 
which is valid also when N = 2. 



Let us see now what problem ( 1.1 ) reads like for radial solutions. Recall that the Riemannian 
Laplacian on the hyperbolic space, for radial functions v = v{r), takes the form 

1 



Av(r) 



(sinh r) 



[(sinhr)^-^ v'{r)]' = v"{r) + {N - l)(coth r)u'(r) , (1.17) 



where the apex ' stands for derivation w.r.t. r. From (1.17) we have that studying energy 



solutions to (1.1 ) for radial initial data is equivalent to studying energy solutions to the problem 

'u^ = (u™)" + (iV- l)(cothr)(u™)' in (0,oo) x (0,r) 
(n'")' = on{0}x(0,T) . 



no > 



:i.i8) 



on [0,oo) X {0} 
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The fact that there exists a finite extinction time T > is a straightforward consequence of 
the vahdity in H^^^iM.^) (indeed also in H^{M'^)) of both a Poincare and a Sobolev inequahty 
(see, for instance, |34l Sect. 5.10] and |7, Sect. 3]), that is 

Ibib < Cp ll^^'lb ) ll'^ll 2JV < Cs \\v'\\2 (1-19) 



iV-2 



for all V G H^^^(M^) and suitable positive constants Cp = Cp{N), Cs = Cs{N), where 

||^;||P= / |7;|P(s) (sinhs)^-Ms, ^^{m^) = {v radial: \\v\\p<oo}. 
Jo 

Moreover, one can prove H Th. 3.1] that the embedding of H^^^{m^) into Lf^j(]H^) is compact 
for all p G {2, 2N/{N — 2)). Notice that, since m G (?7is, 1), this means in particular that 

m + l 

HlM'')<^h^ (H^), (1-20) 
a crucial fact that we shall exploit in the next section. Recall however that the compact 
embedding (1.20) fails in H^{W^), another nontrivial issue that points out the advantage of 



working in the radial framework. 

In the sequel, for notational simplicity, we shall write instead of L'^i^^) and do the same 
for all the functional spaces involved. 

1.2. Plan of the paper. All the above results will be proved in several intermediate steps. 
Local uniform convergence of u™/(l — t/T)™/*^^"™^ to the stationary solution V is shown in 
Section|2| along lines similar to the ones j2]. Then, a suitable upper bound for solutions (holding 
as r — )• oo) is shown in Section [3| whereas a matching lower bound is proved in Section |4] The 
more delicate issue, namely the passage to the relative error u™" /[(l-t/r)"^/(i-"')T/]-l, is dealt 
with in Section [5] Section [6] contains the proofs of the results concerning space-time derivatives 
of solutions, which exploit both regularity theory and the claimed convergence in relative error. 

2. Local uniform convergence of the rescaled solution to the stationary 

PROFILE 



As previously mentioned, each solution to ( 1.1 ) extinguishes in a finite time T > 0. Therefore 



the asymptotic behaviour of u is, from this point of view, trivial: the solution goes to zero as 
t ^ T. In order to study finer properties of u it is very useful to look for separable solutions 



to (1.1) (if any), so that their asymptotic behaviour might unveil at least the expected order 



of convergence to zero of a generic solution. Hence, let us set n(x, t) = g{t)V^/'^{x). After 



some straightforward computations one gets that u is a solution to (1.1) for some uq > (not 
identically zero) if and only if g satisfies 

1 

/ + \ l-m 

g{t)=[l--\ VtG[0,T] (2.1) 



and y is a positive solution to the elliptic problem ( |1.4[ ) for some parameter T > (the 
extinction time). When m G (mg, 1) existence and uniqueness of such a V and its dependence 
the sole radial coordinate r is guaranteed by compactness and by a moving plane method (see the 
fundamental paper [29j). Local regularity and strict positivity of V are instead a consequence 
of standard elliptic arguments. So the velocity of convergence to zero as t \ T for separable 



solutions is given by (2.1). This suggests that, in order to analyze a nontrivial asymptotics, 
it is convenient to study the behaviour of the rescaled solution u{r,t)/g{t). Notice that, if 
u(r,t) = g{t) ^"'^/'"(r), then such rescaled solution trivially coincides with V^/^. For a generic 
u this is of course not true: however, y^/*" seems to naturally maintain the role of an attractor 
for u/ g. 

Motivated by the discussion above, given the extinction time T associated to the solution u 



of (1.1), let us consider the rescaled solution w defined as 



w{r,T)=[^—X'^u{r,t) = e(^~^)Tu(r,T-Te-i), T = T\og(^—^] (2.2) 
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Vr G (0, oo) , \/t£ (0, T) , Vr G (0, oo) . 
Straightforward computations show that w solves the following problem: 



[W 



w 



^(^"') + (T^^^ in (0, oo) X (0, oo) 
on {0} X (0,oo) 



(2.3) 



lio > 



on [0,oo) X {0} 



The aim of this section is to prove that w'^ converges locally uniformly in {r G [0, oo)} to V as 
r — 7- oo (since V is positive, this is equivalent to claiming that w converges locally uniformly 
to y^/™). The basic estimates one needs to exploit in order to prove such result were obtained 
in a celebrated paper [2] by Berryman and Holland, though for regular solutions to the FDE 
on regular bounded domains of M^. Here, first we shall only point out how their techniques, 
with minor modifications, can be applied to this framework too. This will ensure local uniform 
convergence at least away from {r = 0}, while some further work will be required to extend the 
result to neighbourhoods of the origin o. 

To this end, it is convenient to see u as a monotone increasing limit of the sequence of 
solutions {un\ (with extinction times T„) to the problems 



[u 



nit 



Un = 
.mv _ Q 



[u: 



n 



n, 



n 



A«) in (0,n) X (0,r„) 
on {n} X (0, Tn) 
on {0} X (0,r„) 
on [0, n] X {0} 



(2.4) 



non > 



where {non} is a sequence of regular data such that UQnin) = 0, (?io^)'(0) = 0, uon < "^o, which 
suitably approximates uq, and Tn f T. We shall identify n„(-, t) as functions in the whole [0, oo) 
by extending them to be zero outside [0,n]. Notice that (2.4) corresponds to the radial FDE 
with homogeneous Dirichlet boundary conditions posed on the ball of radius n of centred 
at X 



0. 



Lemma 2.1. There exists a positive constant C = C{m,N) such that 

1 

C{T- t)^ < ||u(t)|L+i < 1 - ;^ ll«olL+i Vt G (o,r) 



Moreover, the ratio 



(^x"^)' 



\u 



(2.5) 



(2.6) 



m+l, 



is nonincreasing along the evolution. 



Proof. The left inequality in (2.5) can be proved exactly as in Lemma 1] since it is a direct 



consequence of the Poincare-Sobolev inequalities in (1.19). To justify the right inequality and 



that (2.6) does not increase, it is convenient to pass through the approximating solutions Un- 



Indeed the proof of O Lemma 2] requires the finiteness of the quantity 



/•oo 

/ u"'^\s,t) [A(M™(s,t))]2 (sinhs)^-Ms 
Jo 



for all t G (0, T), which a priori may not hold here. However, one obtains the claimed inequalities 
for Un (we can assume that is regular enough up to the boundary) and then passes to the 
limit as n — )• oo. This is feasible since M^(t) converges weakly in H^^^^ to u"^{t), and by 
monotonicity Un{t) converges to u{t) in L^J^^ and T„ converges to T. □ 



The next result is a key one in order to establish the mentioned convergence of 
stationary profile V. 



to the 
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Lemma 2.2. The following inequality holds true for all r G (0, oo); 



[( 



+ m 



1 



Jo 



m 



;i -m2)r 



w 



m+l 



is,r) 



(sinh s 



ds 



w"'~'{s,a) [wr{s,a)Y (sinh s)^-Ms da 



(2.7) 



< 



m 



m+l 



is) 



(sinh s 



,N-1 



ds. 



Proof. For the smooth rescaled solutions Wn inequahty (2.7) is in fact an equahty, since by 
straightforward computations one verifies that 

d 

d^ 



m 



-W. 



m+l 



(sinhs)^-Ms 



(2.8) 



/ ^(■S)'^) [(^n)r(s,T)]^ X(0,n)('5) (sinhs)^ Ms . 



To get estimate (2.7) it suffices to integrate (2.8) from to r and let n — )• oo: to the first 
integral on the l.h.s. of (2.7) we can apply the weak convergence of u'^{t) to u"^{t) in H^^^ and 
the strong convergence of Un{t) to u{t) in L^^^ , while the second integral is handled by means 
of Fatou's Lemma (thanks to local regularity we can assume that, up to subsequences, {un)t 
converges point wise to ut) or by the fact that 



m + l 
, 2 



m+l 
U 2 



in H\0,T;Ll^). 
Now we are able to prove the following important result. 



□ 



Lemma 2.3. Let w be the rescaled solution (|2.2|) to (|2.3|) and V the radial, positive energy 

(2.9) 



solution to the stationary problem (1.4). Then 



lim ||i«™(r) - V\\^ 



0. 



that is w'^{t) converges uniformly to V in any compact set K d (0, oo) as r —t- oo. 



Proof. We adapt the proof of [2, Th. 2]. First of all notice that, from ( |2.7[ ), one deduces the 
existence of a sequence r„ — )• oo such that 



w"'-'{s,Tn) [Wr{s,Tn)Y (sinh s)^-Ms ^ . 



(2.10) 



This is an immediate consequence of the fact that, thanks to (2.5), the first integral on the l.h.s. 
of (2.7) is bounded from below and the r.h.s. is a constant; therefore the integral 



JO 



w 



m—l 



(s,cr) [wr{s,a)]'' (sinhs)^ ""^dsdo" 



must be finite. Still from (2.7) and (2.5) one gets the boundedness of ||w™(r„) ||2,2; hence, 
up to subsequences, w^ijn) converges weakly in H^^^ to a certain function R. From the 

compact embedding (1.20), such convergence is in fact strong in L^^^"^^^"^ ■ In particular, ii is a 
nonnegative non-identically zero function (indeed (2.5) prevents ||tt;'"(T„)||(m+i)/m from going 

to zero) belonging to -ff^d- 

The next step is to show that R solves (1.4). To this end, take any test function : [0, oo) — )• M 
with compact support in [0,oo), multiply by it the first equation in (2.3) (evaluated at r = t„) 
and integrate by parts in [0, oo). This leads to the identity 



Wris^Tn) 4>{s) [smhs) ds 



{w^)'{s,Tn)(t}'{s) (Sinhs 



,N-1 



ds+ 



(2.11) 



?i'(s, r„) (;/)(s) (sinh s)^ "^ds 



{l-m)T 
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The two integrals on the r.h.s. of (2.11 ) are stable under passage to the limit as n — t- oo: indeed 
w"^{Tn) converges weakly in H^^^^ to R and w{Tn) converges strongly in L^'^^ to B}/'^ (and so 
also in L].^^ locally). Finally, the left hand side goes to zero since its modulus is bounded by 



w 



.1— m 



{s,Tn)(t)\s) (sinhs)^-Ms 



[wr{s,Tn)Y (sinh s)^ ^ds 



which goes to zero thanks to (2.10) and to the boundedness of ||i(;(r„) ||m+i- From the arbitrari- 
ness of (j) and the uniqueness of energy solutions to (1.4) we infer that R must coincide with V. 
Moreover, convergence of w^{Tn) to i? = y in H^^^ is also strong. To prove that, just replace 
(/>(•) by w"^{-,Tn) in the computations above and get 



lim 

n— >oo 



/ [{w^^ns^Tn)]^ (sinhs)^-Ms= / 
Jo Jo 



V "1 (s) (sinh s) 



7V-1 



ds 



(1 - m)T 

[V'{s)f{smhsf-^ds. 

Hence, weak convergence plus convergence of the norms in gives the claimed strong con- 



vergence. Since H^^^ is continuously embedded in (see e.g. Lemma 3.2 below), we have 



proved (2.9) along the special sequence {t„}. To prove that such convergence takes place along 
any other subsequence, one can argue by contradiction. That is, suppose there exists a sequence 
{t/j} such that w'^{tiS) does not converge strongly in to V. By (2.7) we can assume that, 
up to subsequences, w"^{Tk) converges weakly in H^^^ to a certain function Q. Now notice that, 
again from (2.5) and (2.7) (up to a time origin shift), both 



\w{t) 



and 



im+l 
lm+1 



m 



w{t) 



are nonincreasing functions of r. In particular, 

||Q||m+i = lim ||w(rfc)||^ , = lim ||'w(t„ 



im+l 
lm+1 



im 

lm+1 



\V\ 



m + l 



and 



IQ'II < liminf 



m 



[1 - m?)T 



\w{Tk) 



im+l 
lm+1 



(2.12) 



(2.13) 



lim 

n— >oo 



\v'\ 



Since V is the unique minimizer of ||w'||2 among all functions v with prescribed l'^^^'^^^™' norm 
(see |29|). (2.12) and (2.13) necessarily imply that Q = V. Strong convergence of w'^{Tk) to 
Q = V IS then a consequence of (2.13), which leads to a contradiction. □ 



We are left with proving that the uniform convergence (2.9) takes place also down to r = 0. 
In order to do it, we shall use Lemma 2.3 and the following two lemmas, which show how 
positivity and boundedness of w can be extended to a neighbourhood of o. 

Lemma 2.4. For any e > there exist > sufficiently small and > sufficiently large 
such that 

w'^{r,T)>V{Q) - e V(r, r) G [0, r^] x [r^, oo) . (2.14) 

Proof. We can adapt the techniques of |2H Lem. 6.2]. First of all recall that, thanks to the 
local uniform convergence to the stationary profile (2.9), we have uniform boundedness away 
from zero in any compact set which does not contain o. In particular, consider a point xq G 
such that ro = d{xQ, a) G (0, 1/2). For a given tq > 0, set 

k= inf tt;((i(x, o), r) , (2.15) 

(x,r): r>To , x<^B^^/2{xo) 
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where -6^0/2(2^0) is the hyperbohc ball of radius ro/2 centered at xq. Thanks to the observations 
above, A; > provided tq is sufficiently large. Let us consider equation ( |2.3[ ) (more precisely, its 
interpretation as a differential equation on M^) centered at in place of o. To avoid confusion, 
we shall call p the radial coordinate about such xq. Upon defining 

iV=l + (A^-l) sup pcoth(p), 
pe(o,i) 



for any function f = f{p) such that /'(p) < we have 
A/(p) = f"(p) + (iV - 1) cotHp)f{p) > f"{p) + 



iV- 1 



f'ip) VpG(0,l) 



(2.16) 



where the term on the r.h.s. of (2.16) is the Euclidean Laplacian of / associated to the "artificial 



dimension" A^. Therefore in order to seek for a subsolution tp{p,T) to (2.3) centered at it is 
enough to ask (we keep denoting as ' derivative w.r.t. p) 



N - 1 



(2.17) 



as long as p varies in (0, 1). The proof of Lemma 6.2 of [21] ensures that the function 



satisfies (2.17) in the region 



iI>{p,t) = k- 



(P,r)G ( 



{1-p^ 



1 + k 



l-m bp2 



To, To + 



e 

1 — m 



'0 



2 ' y V 4 

upon choosing appropriately the positive parameters /3 = f3{m,N), 6 = 9(171, N) and b = 
b{m, N,k), k being as in (2.15 ). Let us check conditions on the parabolic boundary. For r = tq 



and for p = 1 we have, by construction, ijj = and so trivially ip < w. For p = ro/2 (actually for 
any p G (0, 1)) there holds ip <k, from which ip <why definition of k. Hence, by comparison. 



w {d{x,o),To + rl/A) >V'(d(x,xo),To + ro/4) 



>V(3ro/2,ro + rg/4) = /c- 



ro 



< d{x, xo) < 



(1 + 96A:i-™)i-™ 
3ro 



:Cn>0 



In particular we obtain the existence of a radius ri > and a time ri > such that 

w{r, r) > Co > V(r, r) G [0, n] x {n}. 



(2.18) 



Indeed (2.18) holds for all r greater than ri, rather than only for r = ti. This is a trivial 



consequence of the fact that Co is a subsolution to (2.3) (the comparison condition on the 
lateral boundary {ri} x (n, 00) is satisfied provided Co is small enough, again as a consequence 
of the local uniform convergence (2.9)). 



Finally, we need to refine estimate (2.18). To this end, just observe that the function 



is a solution to the differential equation in (2.3) for any t^: > 0. Still from the local uniform 



convergence (2.9) (and from the fact that V{-) is decreasing) we have that, for any e > 0, we 



can choose r2 = r2{e) < ri and T2 = T"2(e) > ri such that w"^{r2,T) > V{0) — e for all r > T2. 
Therefore g{T), with the choice r.,, = T2, is a subsolution to ( |2.3| ) in the region 

{(r,r) G (0,r2) x {t2,T3)} , 
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where is the time at which g"^{T^) = V{^) — e. Since the constant (V^(0) — e)-^/'" is then a 
subsolution in {(r, r) G (0,r2) x (r3,oo)}, estimate (2.14) follows. □ 



Now we prove the analogue of estimate (2.14) from above. 

Lemma 2.5. For any e > there exist > sufficiently small and > sufficiently large 
such that 

ti;™(r,T) < y(0) + e V(r, r) G [0, r,] X [t„ oo) . (2.19) 



Proof. Again, we shall proceed by constructing a proper supersolution to (2.3). To begin with, 
let a and e two small positive parameters. Our aim is first to obtain a suitable estimate for 
A (V (r/a)) in the region {r < ae}. We have: 



A{V{r/a)) 



1 



,yi/™(r/a) + 



N-l 



'acoth(r) — coth(r/Q)) V'{r/a) 



a^(l — m)T 

The function /i(r) = r coth(r) is regular. In particular, 

h{r) = 1 + h'(0) r + q(r)r^ , h(r/a) = 1 + h'(0)- + qir/a) (-)' 

a \aJ 

where both \q{r)\ and \q{r/a)\ can be bounded by 

h"{s) 



(2.20) 



(2.21) 



Q = max 

se[o,i] 2 



provided a,e are smaller than 1. In order to control the right term on the r.h.s. of (2.20), we 

rcoth(r) — ^coth(r/a) 



use (12.211): 



acoth(r) — coth(r/a) 



h'{0) 



ar 
r 



^'(0) / / X ^ 
+ 2 IV ot)^ 



< 



C 



a a 

for a suitable constant C > independent of a, e. Notice that, since V is regular and l^'(O) = 0, 
there exists D > (independent of a,e) such that |y(r/a)| < De for all r < ae. Hence, 

(iV - l)CDe 



N-l 



a" 



^acoth(r) — coth(r/a)) V'{r/a) 



< 



(2.22) 



Then recall that V{-) is decreasing and V{0) > 0, so that from (2.20) and (2.22) we can claim 
that there exists a constant E > such that for any e > sufficiently small (depending only 
on V, m and A^) there holds 



A{V{r/a)) < 



a 



E 



Vr e (0, ae) 



(2.23) 



From the L'^-L°° smoothing effects (see |21i, Lem. 6.1] or [51 Th. 4.1], together with (2.5)) we 
know that there exist A > and tq > such that u!{r, t) < A for all (r, r) G (0, oo) x (tq, oo). 
Let Ai > ^2 be two given positive constants and let e > be so small that (2.23) holds. For 
a fixed r* > tq + 1 set /(r) = (r — tq)/{t* — tq). First we shall prove that if a > is small 
enough then the function 



ip{r,T) = [Ai{l- f{T)) + A2f{T)]^V{r/a)^ 



is a supersolution to (2.3) in the region 

{(r,r) G (0,ae) x (ro,r,)}. 

To this end, notice that 



(2.24) 



(^(r,ro) = y™(r/Q) > ylf F™(e) VrG(0,ae), (/'(a^, •^) > ^2" l*^) VrG(ro,r*), 

(2.25) 
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while derivatives of ip give 

^r{r,T) = -l/'(r)y^(r/a)(Ai - A2) [Ai{l - /(t)) + A^fir)]^ 



i-i 



> 



m 

^ {Ai-A2)A^~W^{0). 



(2.26) 



m 



A(^-)(r, r) = - fir)) + A2/(r)] A(y(r/a)) < 



A2E 

2j< 



a 



and 



(2.27) 

Collecting (2.26)-(2.27) we get that for to be a supersolution in the region (2.24) it is enough 
to ask 

,:^-l_i.„s . A2E , AfV^{<d) 



1 



(^1-^2) A'" ^-(0)> 



+ 



(2.28) 



m ■ ' " ' a^T (1 — m)T ' 

which is achieved by choosing a = aiV, T, m, N, Ai, A2) sufficiently small. 

Now fix e > sufficiently small. Set Ai = A"^/V{e), A2 = V{0)/V{e), where we assume 
without loss of generality that A™ > V{0) and pick a(y, T, m, N, A, e) that complies with 
( 2.28D . Thanks to ( |2.25[ ) we have 

^{r,To)>A VrG[0,ae], (/?(«£, r) > T/™ (0) VrG[ro,r,]. (2.29) 



From the fact that V{-) is decreasing and from the local uniform convergence (2.9), we can take 
To so large that ■w'^{a£, r) < V{0) for all r > tq (notice that (/? is a supersolution independently 
of To,r* provided r* — tq > 1). Since (2.29) holds we can conclude, by comparison, that w < ip 
in the region (2.24). In particular, 

y2(0) 



V{e) 



Mr G [0, ae] . 



By the remarks above this last result is actually valid for all t^, > tq + 1. Hence, since V{e) — )• 
V{0) > as e — )• 0, we conclude that for any e > there exist so small and so large that 
( |2A9l ) holds true. □ 



Thanks to Lemmas 2.3 2.4 and 2.5 we can extend the result of Lemma 2.3 down to r 
and get the following 







Proposition 2.6. Let w be the rescaled solution (|2.2|) to (|2.3|) and V the radial, positive energy 

(2.30) 



solution to the stationary problem (1.4). Then 



lim Ww^'ir) - V\\l. 



0. 



that is w^{t) converges uniformly to V in any compact set K d [0, 00) as t ^ 00. 

3. Estimates from above 

The goal of this section is to bound the ratio w'^{r, T)/V{r) in L°°((0, 00)) (and not only in 
LJ^^([0,oo)) as we did in Section [2| from above. Since V{r) behaves like e^*^^""'^^'" at infinity 
(see (1.3) or Lemma 3.3 below), it will be enough to give an upper bound for w"^{r, t) / e~^^~^^'^ . 
In fact our goal is to prove the following result. 

Proposition 3.1. Let w be the rescaled solution associated to a nonnegative energy solution u 
to (1.18), as in (2.2). Then, for any e > 0, there exists a positive constant Q' = Q' {uQ,m, N, e) 

N-l , , 

(3.1) 



such that 



w{r, r) < Q' e~^s^ ^ V(r, r) G [0, 00) x [e, 00) 
To this end, we begin with some preliminary lemmas. 
Lemma 3.2. Let v £ H^ad- '^^y ''0 > there exists a positive constant C{ro,N) such that 

v{r) <C{ro,N)\\v\\2,2e Vre[ro,oo). 
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Proof. Consider the function z(r) = {smh.r)^~^ v'^(r). We have: 

z'{r) = (A^ — l)(coshr)(sinhr)^~^ f^(r) + 2(sinhr)^~^ r;(r)w'(r) ; 
integrating between r and tq gives 



z(r) =(sinhro)^~^w^(ro) + (iV - 1) / 7;^(s) (cosh s)(sinh s) 



\Af-2 



ds 



+ 2 / v{s)v'{s){smhs)^-^ds 



(3.2) 



Since ro > and the behaviour at infinity of sinhr and coshr is the same, we can control the 
last two terms on the r.h.s. of (3.2) with a constant (depending on rg and A^) times ||f ||2 2- 
for the first term, notice that H^{ro/2, 3ro/2) is continuously embedded in L°°(ro/2, 3ro/2) and 
H^g^^ is in turn continuously embedded in iJ^(ro/2, 3ro/2) (again, through constants depending 
on ro and N). Hence, there exists C(ro, A^) > such that 

z{r) <C\ro,N)\\v\\l2 Vre[ro,oo), 

which gives the claimed result since sinhr x for r large. □ 

Lemma 3.3. For any fh G (m^, 1) there exists a solution V to 

- AVir) = (r) Vr G (0, oo) (3.3) 
which is smooth, strictly positive, belongs to and satisfies 

A'^ e'^^~^^^ < V{r) < A e-^^-'^> Vr G [0, oo) (3.4) 
for some positive constant A = A{m.). 

Proof. As already mentioned in the Introduction, see |29) and [?]. □ 



Lemma 3.4. Let u be a bounded, nonnegative energy solution to (1.18). There exists a positive 
constant Cq = Co{uQ,m, N) such that 



uir, t) < Co e" ^ V(r, t) G [0, oo) x (0, T) . 
Proof. It is a matter of straightforward computations to show that 



>0 VrG(0,oo) 



(3.5) 



(3.6) 



Thanks to Lemma 



3.2 the fact that ■u™ G -f^^ad ^^"^ boundedness of u, one can choose Cq so 
large that Co e~ 2m is above u on a parabolic boundary of the type [ro, oo) x {0}U{ro} x (0, T), 
for a given ro G (0,oo). The conclusion then follows from (3.6), the comparison principle and 
again the boundedness of u. □ 

We are now ready to prove a better (spatial) estimate from above for nonnegative energy 
solutions to (1.18). 



Lemma 3.5. Let u be a nonnegative energy solution to (1.18). For any G (0,T) there exists 
a positive constant Q = Q{t^:,uo,m, N) such that 

N—l 

u{r,t^) <Qe-^'' VrG[0,oo). (3.7) 

Proof. We shall proceed by constructing a proper barrier. In particular, we shall prove that for 
a suitable choice of the parameter ^ > 0, the following function is a supersolution to (1.18) in 
the parabolic domain (^,oo) x (0, t*): 



n(r,t) = Co Ae-^^V{r)f{t) + e-^''- {I -fit)) 



(3.8) 



where V is the solution to (3.3) associated to a fixed rh G (2m/(l + m),l) (A being the 
corresponding constant that appears in (3.4)) and f{t) : [0,t=i,] — ?■ [0,1] is a regular increasing 
function such that /(O) = 0, /(t*) = 1, which we shall define later. The constant Cq is the one 
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from (3.5): indeed, thanks to smoothing effects (recall the discussion in Section 1.1), there is 
no loss of generality in assuming that mq is bounded. Since 



n(r, 0) = Coe 2m > uo(r) VrG[0,oo) 



and 



uitt) = Co[Ae^^V{0 fit) + e-^^ (1 - /(t))] > Co ^ > u{C,t) VtG [0,t.), 

in order to prove that u(r, t) < u{r, t) for all (r, t) G oo) x [0, t*] we are left with showing 
that, by choosing appropriately ^, there holds 

ut{r, t)>A [u'^ir, t)) V(r, t) G {C, oo) x (0, U) . 



We have: 

ut{r,t) = Cof{t) \Ae^^V{r)-e-^ 



1 r 
m 



Ae—^V{r)f{t) + e-—'- (1 - /(*)) 



m 



--1 



Also, thanks to ( |3.3[ ), 



Ae^^V^{r)f{t) + 



(AT- 1)2 _^ 



e — ^^(2cothr-l) (1 - /(t)) 



1 JV-1a JV-1 



(AT- 1)2 _^ 



e-^'^ (l-/(t)) 



In particular, there exist two positive constants Bo{m,m) and Bi{m,N) such that 



-^2t(r,t) < Co Bo fit) 



■ (JV-l)(l-m) ^ (iV-l)(2-m) „ „ J_ 

e 2m e 



L)(2 — m) 1 -, A^— 1 1 1 

— f"\t) + e-^'' (1 - f{t))^~^ 



-A(n™(r,t)) >Co™i?i 



JV-l f JV-1 , , JV-1 

e 2 « e 



f{t) + e-^^ (l-/(i)) 



Hence it is enough to show that, if ^ is properly chosen, the following inequality holds in 
(^,oo) X {0,U): 



(JV-l)(l-m) c (JV-l)(2-m) „ „ 



e 2m ^ e 

JV-1 



'i^^'-f^-\t) + e"^'' (1 - /(t))- 



(3.9) 



Ccl-'"i?*/'(t) 

< 6-^^7(0 + e-^'- (1 - fit)) , 

where 5=,, = B^{m, fh, N) is another suitable positive constant. Applying the change of variable 
p = r — ^ and using the fact that m < 1 we infer that (3.9) is implied by 

(JV-l)(2-m) „ 1 



Co^-"i?./'(t)e-^« 

(JV-l)(2-m) 



-it — -PfT^~^{t) + e-^P (1 - /(t))' 



-1 



< e 2m 



f{t) + e-^P (1 - /(O) V(/5, t) G (0, oo) X (0, t.) 



(3.10) 



Since we choose m to lie in the interval (2m/ (1 + m), 1) we have that 

e = {N-l) 



therefore (3.10) reads 



Cl-'^B, e-^« f'(t) 



1 1 1 ^ 

2m fh 2 



-'^^'"^^^-^(t) + e-'^P (1 - /(t))^-^ 
e-^P/(t) + e-^7l-/(t)) 



< 1. 



(3.11) 



Our aim is now to show that for a suitable choice of f[t) the function /(p, t) remains bounded in 
(0, oo) X (0, t^). To this end, let us set x = e~^'^~^^P and f(t) = h{t/t^,), h being a function to be 



14 



GABRIELE GRILLO, MATTEO MURATORI 



defined which satisfies all the fulfilments required to / in the interval [0, 1] instead of [0, t^,]. The 
boundedness of l{p,t) is implied by the boundedness of the ratio (recall that m > 2m/(l + m)) 



i'^) / 1 .i\ , . —. — = (^)- 



xii^+^)h{T) + x5 (1 - /i(t)) 



x^/i(r) + 1 - /i(t) 



(3.12) 



a 



m 



for (x,r) G (0, 1) x (0, 1). If a > 1, which corresponds to m < 1/2, the numerator in (3.12) is 
always smaller than or equal to the denominator. Therefore we remain with the case a G (0, 1), 
that is m E (1/2, 1). Here it is convenient to choose h as 



Kr) = 1 - (1 - r) 



1 

T I « 



In this way, performing the change of variable a = 1 — t, (3.12) becomes 



— (J a 

" X 2 ( 1 



1 

(J a 



+ X 2 (J 



(3.13) 



(T OL 1 + CJ a 



for {x,a) G (0,1) x (0,1). Now notice that, for a varying in [1/2,1), (3.13) is bounded by a 
constant that depends only on a (recall that a G (0,1)). If instead a varies in (0,1/2), the 
boundedness of (3.13) is equivalent to the boundedness of 

-1 



X a + x2 

X 2 cr a + 1 



< 



X a 



a + l _ 1^ 
X 2 a a+l 



+ 1. 



(3.14) 



For any fixed x E (0,1), the maximum of the r.h.s. of (3.14) as o" G (0, oo) can be found 
explicitly, and it is equal to 

xt(i-°)a"(l-a)^-" + l. 

Summing up, we have proved that for a suitable choice of f{t) (depending on whether m < 1/2 
or m G (1/2,1)) the function l{p,t) in (3.11) is bounded by a positive constant K = K(tif,m) 
as {p,t) varies in (0,oo) x (0, t*). This means that if we set 

e> - log {C^-^B^K) 



we ensure that u{r,t) as in (3.8) is a supersolution to (1.18) in the region (C,c«) x (0,t*). By 
the comparison principle, in particular. 



u{r,U) < u{r,U) < Cq A^^ e 2" 



Vr G (^, oo) 



Since u{-,t.t) is also bounded in (0,^], this gives (3.7). 



□ 



The result just proved is not enough in order to bound from above the ratio w"^{r, T)/V{r) 
since it provides such boundedness only at any fixed r G (0,oo): indeed, recall that u{-,t) is 
bounded (so is 'w{-,t)), V{r) is locally bounded away from zero and its behaviour at infinity is 
the same as e~^^~^^^'. What estimate (3.7) lacks is a decay rate of order (T — t*)^/^^"'") on the 
r.h.s.: we shall now prove that this is indeed the case. 

Lemma 3.6. Let w be the rescaled solution associated to a bounded, nonnegative energy solution 
u to (1.18). There exists a constant Ci = Ci{uo,m, N) such that 

N~l 

w{r,T) <Cie-^'' V(r, r) G [0, oo) X (0, oo) . (3.15) 

Proof. The fact that 

\\um^<D{T-t)T^^ VtG(0,r) (3.16) 

for a suitable positive constant D = D{uQ,m, N) can be proved exactly as in |2H Sects. 5-6] 
(see also the results of [8J). Then notice that Lemma 2.1 in particular yields 

||(^")'(r)|L< 



Vr G (0,oo) . 



(3.17) 
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The Poincare inequality in (1.19), Lemma 3.2 and (3.17) provide the following bound: 

N-l 



w 



\r,T)<D'{ro,N) \Uu^)'^ 



V(r,r) G [ro,oo) x (0,oo) , 



which together with (3.16) gives the claimed estimate (3.15). 



□ 



Notice that we proved (3.15) under the hypothesis that u is a bounded, nonnegative energy 
solution. However, thanks to the aforementioned smoothing effects, it also holds for any solution 
associated to data as in the hypothesis of Theorem |1.1[ provided one starts from t = e > 
instead of r = 0. 

The bound (3.15) for w is the exactly the same as ( |3.5| ) for u, which was a key starting point 
in order to prove the claim of Lemma 3.5 Indeed, as we shall see now, the barrier exploited in 
the proof of such Lemma also works for the equation solved by w. This allows us to obtain the 
next fundamental estimate from above. 



Proof of Proposition 3.1 As just remarked, we need only prove that for a suitable choice of the 
parameter ^ > the barrier constructed in the proof of Lemma 3.5 still works. So, for a given 
G (e,oo), consider again the function 



w{r, t) = Ci a e V? V{r) fir) + e" (1 - /(r)) 



(3.18) 



where ^, V, A, f are as in (3.8) (just replace f by r and let /(e) = 0) and Ci is the constant 
appearing in ( |3.15 ). First of all, since /(•) is always included in [0,1], m < 1 and m G 
(2m/(l + m), 1), we can bound the reaction term in (2.3) (what makes (2.3) actually differ from 
( |1.1[ )) in the following way: 



1 



(1 -m)r 



w{p + £„t) < CiB2e~ 



N-l „ 1 , , iV-1 „ 



_ iV-1 



7-(T) + e- 
7(r) + e- 



(l-/(r))' 
l-/(r))" , 



where B2 is a positive constant depending only on T, m, m and we performed the usual change 
of space variable p = r — ^. So, the equivalent of (3.10) in this context reads 

(JV-l)(2-m) 



Cii3o/'(T)e-^ 



^ ^/™' (r) + e-^P (1 - /(r))' 



f{r) + e-^ 

V(/),t) G (0,00) X {e,n) 
By elementary computations one gets that 



(1 - fir)) 



- < CrBi-CiB2e-'^ 



provided 



e > - log 

e 



2Cl-"'B2 

b[ 



(3.19) 



Therefore, under assumption (3.19), we can repeat the same proof of Lemma 3.5 starting from 
(3.11) (one replaces B^, by 2Bq/Bi and Cq by Ci). Hence we end up with the existence of a 
positive parameter = S,{Ci{uq, m, A^, e), T, e) such that 

w{r, t) < w{r, t) V(r, r) G [C, 00) x [e, r*] . (3.20) 

The validity of ( 3. 1 ) is then a consequence of (3.20 ) , ( 3.18| ) (evaluated at r = t^,) and ( 3. 16 ) . □ 



Remark 3.7. In the proofs of Lemmas |3.4[ |3.5| and Proposition |3. 1| we applied the comparison 
principle in parabolic regions of the form (,^, 00) x (0, t^,), without considering the {r = 00} side 
of the parabolic boundary. However, this technical issue is easy solvable by approximating u 
with the solutions u„ to (2.4), applying comparison between tt„ and u in (^,n) x (0,t*) and 
passing to the limit as n — )• 00 (using also that T„ f ^)- 
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Remark 3.8. The estimate from above provided by Lemma 3.5 actually holds for all m £ (0, 1), 
since its method of proof only requires the existence of a solution to (3.3) satisfying (3.4) for a 
value of fh which can be taken as close to 1 as necessary. 

4. Estimates from below 

This section is the dual of the previous one: our aim here is to bound the ratio w'^{r, r) /V{r) 
from below. Again, thanks to (3.4), this will be equivalent to providing a lower bound for 
w"^ {r, t) / e~^^~^^ ^ . More precisely, we shall prove the following result. 



Proposition 4.1. Let u be a nonnegative non-zero energy solution to (1.18). For any e E (0, T) 
there exists a positive constant P' = P'{uo,m,N,e) such that 



N~l 1 

uir,t^) > P' e-^" {T -t^)^ 'i{r,t^) € [0,oo) x [e, T) . 



(4.1) 



To this end, we need a preliminary step. 



Lemma 4.2. Let u he a nonnegative non-zero energy solution to (1.18). For any given a > 
N — 1 and t-t G (0, T) there exists a positive constant P = P{tif, a,UQ,m, N) such that 

u{r,t^) > Pe-^'' VrG[0,oo). (4.2) 

Proof. We shall prove that for a suitable choice of the positive parameters fiQ and ^ the following 
function is a subsolution to (1.18) in the parabolic region (^, oo) x {t^,/2,t^,): 



n(r,t)=/.o[(l + e-"(^-«))/(t)-l 



(4.3) 



where / : [t^/2,t^] — )• [1/2,1] is an increasing function such that f{t^,/'2) = 1/2, /(t*) = 1, to 
be defined later. We have: 

1 



n,(r,t)=/xo-/'(t)(l + e-"(^'-«) 



m 



1 + e 



fit) - 1 



and 



A «(r, t)) = fil^ fit) e-"('^-«) (q2 - (AT _ l) acothr) 

in the region where (l + e"" (^-«)) /(t) - 1 is nonnegative (below, we shall always work tacitly 
in such region), while both and A(«™') are zero outside it. Let us check conditions on the 
parabolic boundary. On [^,oo) x {t^:/2} u satisfies 



1 



1 + e 



{r-0 



1 







u{r,t^/2) = /zo 

and on {S^} x (t^,/2,t^,) there holds 

u{C,t) = fio [2f{t)-l]^ <^o. 

Therefore, in order to have u < -u on [^,cxd) x {t*/2} U {^} x (t,,/2,t*), we need only ask 

fio < inf n(^,i) = A(^,t*,tio) > 0, 
ie(^,t.) 

the last inequality following from standard positivity results (if t^, is close enough to T, one 
can also exploit the results of Section |2]). Now let us check the differential equation. Upon the 
usual change of spatial variable p = r — ^, 



Utir,t)<A{yrir,t)) V(r, t) E (^, oo) x 



1 1^ 



reads 



1 



1^0 -fit) (1 + e-"'') [(1 + 
<MS^/(t)e-"'' 



-a p 



) m - 1] 



- (A^ - 1) a coth(p + 0) V(p, t) G (0, oo) x ( ^ , 



(4.4) 
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If we choose ^ = S^{a,N) so large that, for instance, 



1 



a 



(TV - l)acothC > - (a^ - (N - 1) a) = Ci{a,N) > 0, 



(4.5) 



we get that (4.4) is imphed by 

(l + e-"0 [(1 + e-" ")/(*) -1]' 



1 



fit) 



mCi ' ' f{t)e-'^P 
which is in turn imphed by (recall that 1/2 < f{t) < 1) 



< 1 V(p,t) G (0,oo) X ( ) , 



m Ci 



f'{t)e-P [(1 + e- 



-^)/(t)-l]^-'<l V(p,t)G(0,oo)x (4.6) 



L{p,t) 



If m < 1/2 the function L in (4.6) is bounded from above (take / regular) by a constant that 
depends only on t^,, a, m and A^. If instead m S (1/2, 1) this is in general false, unless one 
chooses / carefully. To this end, consider the function 

1 



Mr) 



Vr G [0, 1] 



and set 



fit) = h{2--l 



Elementary computations (one can find the exact maximum of L{p,t), in the region where 
(1 + e~°'P)f{t) > 1, at any given t) show that 



L{p,t)<— V(p, t) G (0, oo) X 



" ; t^ 



for a suitable positive constant C2 = C2{a,m, N) (which we assume to work for the case 
m < 1/2 too). Hence, we proved that u as in (4.3) is indeed a subsolution to (1.18) providing 
that 



fio < min A(^(q;, N),U,uo) , 



C~2 



1 > 

1 — m 



In particular, at t = there holds 

u{r, t^) > u{r, t*) = po (t*, a, uq, m, N) e™ e~ ™ Vr G (^, 00) , 

which yields the thesis together with the local positivity of ti(-, t,,,) in (0, ^). □ 

The result provided by the previous Lemma asserts that at any given t^: G (0, r(no)) all 
nonnegative non-zero energy solutions of (1.18) go to infinity (as r — ?■ 00) slower than e"™*", 
for any a > N — 1. Exploiting such fact, we are able to prove that one can actually take 
a = N-l. 



Proof of Proposition 4 ■ 1 Once again we construct a lower barrier which has the desired property 
stated in (4.1). Indeed, given G (2e,T) (the final result will follow just by replacing e with 
e/2), consider the function 



u{r, t) = po 



e-Hr-O + A-V(r)e(^-i)«) f{t) - e'^^^'^) 



+ e 



where a = a{m, N) and f3 = f3{m, N) are fixed parameters such that 
a> N -1, (3 <N -1, a<l3+{N-l) 



m 



(4.7) 



(4i 



/ : — >• [0, 1] is a regular increasing function that satisfies /(e) = 0, /(t*) = 1 and V is the 
solution to (3.3) associated to fh = m {A being the corresponding constant appearing in (3.4)). 
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We want to prove that, if one chooses the positive parameters ^ and fiQ properly, then u is a 
subsolution to (1.18) in the parabohc region {£,,00) x (e,t*). By (4.7), we have 

u{r, to) = /io ^ e~ ™ Vr G 00) 

and 

Therefore u and u are correctly ordered on the parabolic boundary [,^, 00) x {e} U x (e, t^<) 
provided (recall (|42|)) 



l_iQ < e P{e,a,uo,m, N) and /Uq < 2 ™ inf u{£,t) = X{£,t^:,uo,m, N,e) > . 



that is 

/io < min ^e^ ™ ^ P , A 
Now let us compute the derivatives of u{r,t): 

u,{r, t) =/xo - fit) (e-^ + V(r)e(^-i) sign+ g(r, t) 
A«(r,t)) 



(4.9) 



+ 6" 



(/(t) - 1) (/32 - (iV - 1) /3cothr) e-^(^-«) - ^'Vsj(r)e(^-i) «/(*)] sign+ g(r, t) 
+ (a^ - (AT - 1) a coth r) e"" , 



where for the sake of notational convenience we have set 



If we take ^ = ^(a, A^) large enough so that ( [4.5^ holds, set p = r — £ and use (3.4), we obtain: 
1 • ^ 



UtiP + £,t)<i,o- fit) (e-'3^ + e-(^-i)^) [(e 



)/(t)-e-^^' 



+ e 



-a p 



A(n"(/9 + ^,t)) >/x[r [-A^e-(^-i)^«/(t)e-^^ + Cie-"'' 

V(/>,t) G (0,00) X (e,t^) . 
Hence once we choose = ^(a, m, A^) such that, in addition to (4.5 ), it also satisfies 

- 2 ' 

we deduce that in order to have u < u in (C,cxd) x (e,t*) it is enough to ask (recall that 

/We [0,1]) 



m \ 



e-(iV-i)p + e-"P 



^ PO o ^ 



(4.10) 



y{p,t) G (0,00) X {e,Q. 

Upon setting f{t) = h{{t — s)/{t^. — e)), h : [0,1] — ?• [0,1] being a given regular, increasing 
function such that /i(0) = and h{l) = 1, (4.10) reads 

2 / t-e \ / . .^ . \ r 1:^-1 



h! 



mCi {U — e) \t* — e 



< 1 (4.11) 



L{p,t) 

y{p,t) G (0,00) X {e,U) 
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Thanks to ( |4.8| ) and to the fact that G (2e,T), the function L{p,t) in ( |4.11 ) is bounded in 
(0, cxd) X (e,t*) by a positive constant C2 = C2{T,m,N,e). Therefore u{r,t) is a subsolution to 
( |1.18[ ) providing that (recall ( |4.9[ )) 



A*o < mill e 



inf u{^,t) 

ie(e,i*) 



which is implied by 



fJ'O ^ C's mill ( 1 ) inf t) 

te(£,t,) 



(4.12) 



for a suitable positive constant C3 = C'^{S,,UQ,m,N ,e). Exploiting Proposition 2.6 we can give 



a quantitative lower bound for the r.h.s. of (4.12). Indeed (2.30) yields (in particular) the 
existence, for any given > 0, of a time t = t{^, uq) such that 



u{C,t)>^^iT-t) 



1 

1 — m 



e {i, T) 



(4.13) 



From standard positivity results we also know that, should e < i, u{(^, t) is still positive between 
e and t; this fact and (4.13) (together with the local positivity of V) ensure the existence of a 

\fte{e,T), 

(4.14) 



positive constant C4 = C4{S,,Uo,m, N,e) such that 

u{tt) >C4{T-t)^- 

which gives 



inf u{^,t) >Ca{T-U)t-- 



Combining (4.14) and (4.12) we infer that (4.12) is implied by 

/uo <C5 (r-t,)T^ (4.15) 

for another positive constant C5 = C^{^^ uq, m, N, e). The validity of (4.1 ) for r varying in (^, 00) 
is then a consequence of choosing fiQ as the r.h.s. of (4.15 ) and evaluating the subsolution u(r, t) 
at t = t^:, recalling (3.4). On the contrary, the validity of (4.1) as r varies in (0,^] is a direct 
consequence of (2.30) and the local positivity of u. □ 



Proposition |4. 1| can of course be reformulated as follows. 

Corollary 4.3. Let w be the rescaled solution associated to a positive energy solution u to 
(1.18). For any tq > there exists a positive constant P" = {TQ,UQ,m, N) such that 

V(r,T) G [0,00) X [to, 00). 



wir, t) > P" e 



Notice that, as already anticipated in the Introduction, from Propositions 3.1 and 4.1 



deduce formula (1.9) in Theorem 1.3 



we 



Remark 4.4. Both in the proofs of Lemma 4.2 and Proposition 4.1, when computing the 
hyperbolic Laplacian of the barriers (4.3) and (4.7), we neglected Dirac terms coming out from 
the second derivatives of positive parts. However this makes no problem since it is easy to check 
that such terms are nonnegative. 

Remark 4.5. When we applied the comparison principle to u and u, both in the proofs of 
Lemma 4.2 and Proposition 4.1 , we did not take into account {r = 00} as part of the parabolic 
boundary. To justify more rigorously those passages, it is enough to consider the following 
family of modified barriers: 



e]T <u, e > 0. 



where u is either (4.3) or (4.7). Straightforward computations show that is a subsolution 
to ( 1.18| ) as long as u is. Moreover, for any fixed e > 0, u^{-,t) has the property of being zero 
outside a compact set of the form R{e)] C [£,, 00). One then applies the comparison principle 
in (^, i?(e)) X (e, t*) (let e = t*/2 when u is as in (4.3)), gets u^<u and lets e — )• 0. 
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Remark 4.6. The estimate from below ( |4.1[ ) holds for all m G (0, 1), even though one has to 
admit that the constant P' in there, for m G (0, m^], depends on t* as well. Indeed, in the proof 
of Proposition 4.1 for simplicity we exploited the existence of a solution V to (3.3) satisfying 
(3.4) for m = m. However, if m G (0,ms], such a solution does not exist. Nonetheless it is easy 
to see that the choice of any m G {rris, 1) instead of m would have worked the same (provided 
one requires in addition that a < (N — l)/m). 

Hence, this result together with the one discussed in Remark 3.8 prove the bounds ( 1.13[ ) of 
Theorem Ol 



5. Results for the relative error 



As previously discussed, our interest is to show that the solution u to (1.18) converges in a 
strong sense to the stationary state V which solves (1.4). More precisely, in this Section we 

(5.1) 



shall prove uniform convergence in relative error, namely 



lim 

T— >00 



V 



1 



0, 



that is (1.5) rewritten in rescaled variables. 

Si nce V is strictly positive in any compact subset of H^, as a direct consequence of Proposi- 
tion 2.6 (formula (2.30)) we already know that (5.1) holds locally on . The nontrivial point 



is to prove that (2.30) holds up to r = oo, what we are concerned with in this section. 

: uo"^/V — 1. Using 



To begin with, let us write the equation solved by the relative error 

V'6' 



the fact that V satisfies (1.4), we get: 
1 



(1 + 



1 



(1 + 



[1 + 



(5.2) 



m ■ ■ yti (1 - rn)T L 

where, as remarked above, the apex ' stands for derivation w.r.t. r. It is also important to 
consider the equation solved by minus the relative error ip = 1 — w'^/V, which is 

— (l-V)-"Vr = ^^^^AV' + 2^^+ ^ 

m 



:i 



Vr 



{l-m)T 



(i-V')-(i-^)' 



In order to prove (5.1) we shall construct suitable upper barriers both for (p and ip, following 
the approach of |lOj . 

The next lemma shows a good approximation property for V, which will turn out to be very 
useful to overcome some technical difficulties related to the upper barrier for ip. 

Lemma 5.1. There exists a sequence of positive, regular radial solutions {Vn} to 



-AVn 



1 T/"i 

(l-m)T„ 



(5.3) 



on Bn 

[Vn = on dBn 

where Bn is the hyperbolic ball of radius n centered at a, such that: 

• yn{f) ^ V{r) for all r > 1 and Vn ^ V pointwise; 

• Tn^ T; 

• for e > arbitrarily small, one can choose and r^ so large that the following inequality 
holds: 

V^{r) < {1 - e)V'{r) Vr G {r„n) , > . (5.4) 

Proof. For any given n G N, let Xn be the set of all functions v G -ff^ad such that v{r) =0 for 
all r > n and ||w||m+i = ||y || m+i . Consider the following minimization problem: 



Find Vn G X„: 



n Il2 



mm \\v 

v£Xn 



(5.5) 



Thanks to standard arguments (in particular, the compact embedding of H^^^ into ^rad^^^^™)' 
the solution Vn of (5.5) exists, is unique and solves (5.3) for some T„ > 0. Since X„ C Xn+i, 
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the numerical sequence 



-'nil 2 



} is non-increasing. This imphes that {T„} is non-decreasing: 



indeed multiplying (5.3) by f„ (upon replacing Vn by Vn) and integrating by parts yield 



m + l 

nil"' 
"nil m+l 



\v\ 



m + l 
I m 
I m + l 



"i)IKIl2 (i-"i)IKIl2 



Let us call T the limiting value of {T„}. Because {||i'^||2} is bounded, {fn} (along subsequences) 
converges weakly in H^^^^ and therefore strongly in L^^^^^^"^ to a certain function V. Passing 
to the limit in the (weak formulation of the) equation solved by Vn, we get: 

1 



AV 



(1 -m)r 



V 



"i on 



(5.6) 



First of all note that T cannot be oo: if it were, from the Poincare inequality in (1.19) V would 
be zero, while we know that 



l^ll(m+l)/rn. — 1™ ll'^n || (m+l)/m 



{m+l) /m 



> 0. 



So y is a positive, energy solution to (5.6) for some T < oo, having the same L 



- (m+l)/m 



rad 



norm 



as V: by the uniqueness result recalled in the beginning of Section [2] it necessarily coincides 
with V (and so r„ f T). 

The next step is to prove that {fn} converges to V also in L^[^- Since {vn} converges 
pointwise to V (up to subsequences), and locally L^^^^^^'^ is continuously embedded in L^[^, 
we need only dominate {vn} outside a compact set with a function that belongs to L^^^. To 
this end note that, from Lemma 3.2 and from the fact that {||f^||2} is non-increasing, there 
exists a constant K > such that 



Vn{r)<Ke 2 Vr e [l,oo) , Vn. 



(5.7) 



to L^J^ ■ Therefore we have to improve (5.7) using the equation solved by Vn'- integrating it 



When m is greater than or equal to 1/2, the function on the r.h.s. of (5.7) does not belong 

to L^J^ ■ Therefore we have to improve 
between and r we obtain the equality 

'""vk^s) 



(sinhr)^ ^ v'n{r) 



(sinh s)^-Ms 



(5.8) 



Suppose now that (5.7) holds with a generic exponent —{N — l)a (let a > 0, a 7^ m) in 
place of — (A'^ — l)/2. Exploiting the corresponding analogues of (5.7) and (5.8), after some 
straightforward computations we get: 

_^;(^)<Q^I^e-(^-i)-fe(^-i)^'-_l) VrG[l,oo), (5.9) 



m — a \ 

where Q > is a suitable positive constant that does not depend on n (which may change from 
line to line). Since Vn{n) = 0, we have: 

v^{r)= -v'^{s)ds VrG[l,n). (5.10) 



If a > m (5.9) gives —v'^ (r) < Qe-(^-l)^ in which case an integration of (|5.10|) entails 

Vn{r) <Qe-^^-^^'' VrG [1,00). 



(5.11) 



If instead a < m (5.9) gives —v'^{r) <Qe '■^ and still integrating (5.10) one gets 

Vn(r) < Qe~(^~^)^'' VrG [1,00) 

(the case a = m can be dealt with similarly). It is plain that starting from a = 1/2 and 
proceeding in this way, after a finite number of steps we obtain (5.11). Since e"*-^""'^^'' belongs 
to iJad^' have our dominating function and the convergence of {vn} to V takes place also 
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in L^[^- Such convergence is crucial because it ensures that (5.4) holds for the sequence {vn}- 
Indeed applying (5.8) to Vn = V yields 

lim -(sinhr)^-iy'(r) = H (sinhs)^-Ms = > ; (5.12) 

r->oo Jq T 

for a given e > 0, take so large that 

1 



J '^(sinh.)^-M.<C7y|. 



Thanks to the just proved convergence of {vn} to V in L^[^ and to the fact that T„ f T, there 
exists n = n(r^, e), sufficiently large, such that 



1 



-(sinhr)^-^<(r)= F %^ {smhs)''-'ds > T (sinh s)^-Ms 

Vr G (r^, n) , Vn > n{r^, e) . 

Hence there holds 

-v'Jr) -(sinhr)^-i<(r) Cy(l - e) ^ , . ^ , . 



that is (5.4) for the sequence {vn}- 

In general it is not guaranteed that Vn{r) < V{r) in [l,oo). Therefore, in order to conclude 
our proof, it is necessary to modify {vn}- To this aim, consider the following sequence: 

Vn — An Vn , 

where, for any fixed n. An is the largest number belonging to (0, 1] for which An Vn{r) < V{r) in 
[l,oo). We can assume that An < 1 eventually, otherwise there is nothing to prove since along 
a subsequence {Vn} has all the properties claimed in the statement of the Lemma. If instead 
An is strictly lower than 1, recall that each Vn is compactly supported and V is strictly positive, 
therefore Vn{r) necessarily touches V{r) at some point r = ^n S [l^co) (see figures [l]and[2]). 
Now there are two possibilities: either the sequence {^n} remains bounded or it is unbounded. 
In the first case, along subsequences {^n} converges to a certain value ^ G [l,oo). Since {vn} 
also converges locally uniformly in [l,oo) to V (see Remark 5.2 below), then 

Vn{U^V{0>0; 

but by definition of (^n, 

An^n(en) = nen)^ne). 

Hence An — ?• 1. In the case S^n ^ oo (again, along subsequences), clearly each lies in the 
interior of [l,n) eventually. Therefore, in addition to Vn{S,n) = ^(Cn)i we also have V^('^n) = 

^'(^n). So, 

-An (sinhCn)^-^ vU^n) = "(sinh ^n)"^^' V:,iU) = -(sinh^n)^"! V {(n) ^ Cy > , 

but also 



1 



-isinHnf-'vUCn)= /"^(sinhs)^-M.^ r ^^(sinh.)^-M. = C7y, 

where we have used once again the fact that {vn} converges to V in L^[^ and Tn t T. This 
means that in this case as well {An} is forced to go to 1 and so {Vn} is indeed a sequence which 
has all the required properties (no subsequence is needed since the ongoing proof actually holds 
along any subsequence). Just observe that the parameter Tn appearing in (5.3) here is actually 

r„Ay"^-\ □ 



RADIAL FAST DIFFUSION ON THE HYPERBOLIC SPACE 



23 




Figure 1. Vn touches y at ^„ = 1. 




Figure 2. K touches y at ^„ > 1. 



Remark 5.2. As anticipated in the proof above, convergence of the sequence {vn} solving (5.5) 
to V also occurs locally uniformly in [l,oo). This fact can be proved in several ways. One of 
them is the following: {vn} converges weakly in H^^^ to V (so strongly in L^^^^^^^"^), but using 
the equation solved by each Vn and the fact that T„ f T one gets ||f^||2 II ^'11 2- Therefore 
such convergence is actually strong. Since H^^^ is locally (in (0,oo)) continuously embedded in 
L°°, the assertion holds. 



For reasons that will become clearer later (see the proof of Proposition 5.4 below), instead 
of analysing minus the relative error ip it is convenient to study its natural approximation 
tpn = ^ — w'^/Vn, where {Vn} is the sequence coming from Lemma 5.1 It is straightforward to 
check that ip > tpn and the equation solved by ipn is the following: 



— (l-V'n)" 

m 



i-i V'lb' 1 

Vn ^/^i^n + 2-^ + 
Vn-" 



1 

Tn. 



(1 - ^n) 



(5.13) 



valid in the region {(r, r) E (0,n) x (0,oo)}. 

The following fundamental Lemma provides a family of supersolutions to (5.13). 

Lemma 5.3. Let tq > 0. Consider the function 

^{r,T) = C---A{T-To). 
r 

If the positive parameters A, B, C, f satisfy the condition 

5^ .(^-i)(i-e)(^~i)f 



A 

m (1 



m)T 



(5.14) 



(5.15) 
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where e > is a small fixed number (e = 1/5 will work) and K^^ is a positive constant, then ^ 
is a supersolution to (5.13) in the region {(r, r) : ^(r, r) G [0, 1)} n {(r, r) G (f,n) x (ro,oo)} 
independently of n provided 

f > f V r* , n > n* , 
r* and n* being positive numbers depending only on e, m and N. 



Proof. Let us compute the derivatives of 

B 

r 



-A, ^' 
-BA 



B 



(N-l) 



^, coth(r) 



In order to ensure that ^' is a supersolution to (5.13) it is enough to ask 

1 



m 



that is 
A 



m 



(1 _ > 



B 

T 

T/m 



Vn 



{N-l) 



coth(r) 



1 — m)T„ 



1-^) , 



v 



+ 



1 



(1 - m)T, 



(1-^), (5.16) 



n 



where we have dropped the nonpositive term —(1 — Now we need to suitably estimate 

from above the r.h.s. of ( |5.16 ). To this end, take e > sufficiently small. Moreover, take uq 
and ro so large that T„ > T/2 for all n > uq and coth(r) < l + e for all r > rg. Thanks to (5.4) 
and to the fact that Vn <V in [l,oo), we have: 



TZ < 



B 

T 

T/m 



(iV-l)(l + 6)§ + 25 



+ 



{l-m)T 



<-^[(^-l)(l + e)^ + 2(l-.r]+^(l-^) 

' Vn 

V(r, r) G (ro V V 1, n) x (tq, oo) , Vn > no V . 



Recall that the stationary solution V satisfies (5.12); in particular, by (1.3), there exists ri 
sufficiently large such that 



V'{r) <-{N -l){l-e)V{r) Vr > ri 



Therefore, 



7^<-(iV-l)^[-(l + 6) + 2(l-6)2]+ 2 



m)T 



(5.17) 



V(r, r) G (ro V ri V r^ V 1, n) x (ro, oo) , \/n > uq V n,: . 

Since e is small (it is enough that — (1 + e) + 2(1 — e)^ > 0) and Vn <V in [l,oo), we can replace 
by y in ( 5.17| ), which yields 



n < -^ .(iv-i)[-(i + ^) + 2(i-.n +,^ (1 - *) 



(5.18) 



V(r, r) G (ro V ri V r^ V 1, n) x (ro, oo) , \/n > nQ\/ 



Collecting (5.16) and ( |5.18 ), we deduce that a necessary condition for 4* to be a supersolution 
to (|5.13|) in the region {(r, r) : ^(r, r) G [0, 1)} Pi {(r, r) G (ro V ri V r^ V 1, n) x (ro, oo)} is 



— > ^ (1 - M/)^-^ + ^ (1 - *) 



(5.19) 
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Still from (5.12) and (1.3) we get that 



Vr>r2 



for another positive constant K^^ as long as r2 is large enough. So, the final condition on f is 



m (1 — m)l 



(5.20) 



Summing up, if one fixes the positive parameters e (small enough, say e = 1/5), ^, i?, C, chooses 
r* = ro V ri V r2 V V 1 and takes f so large that (5.20) holds, then the function ^ as in ( 5.14[ ) 
is a supersolution to (5.13) in the region {(r, r) : ^'(r, r) G [0, 1)} n {(r, r) S (f, n) x (ro,oo)} 
for all f > f V r* and n > ng V = n* . □ 

In the next theorem we shall see, along the lines of the proof of [lOj, Th. 2.1], how to use 
the (family of) barriers provided by Lemma 5.3 in order to prove that ^ becomes smaller and 
smaller as r — )• oo. To this aim it is crucial to recall that Theorem 1.3 which we proved in 
Sections |3]-|4] implies the existence of a time > and two positive constants cq, ci such that 

ii;™(r, r) 



V{r) 

Proposition 5.4. Let ij) = 1 



< ci V(r, r) G [0, oo) x [tw, oo) , cq < 1 < ci . 

w'^ /V . There holds 

limsup sup '4){r^T)<{). 

T-s>oo re[0,oo) 



(5.21) 



Proof. First of all consider the barrier ^ as in ( 5.14[ ), with the choices A = B = 1. Then set 
C = 1 — Co/2 and assume that f is greater than (2/co) V f V r*, where f is taken large enough 
so as to satisfy 

m (1 — m)l 

that is (5.15) when A = B = 1. By construction, such a ^' is always lower than 1. Therefore, 
from Lemma 5.3 it is a supersolution to (5.13) in the region (f,n) x (ro,oo) (let tq > t^) for 
all n > n* as long as it is greater than or equal to zero. Since we want to compare ^ with the 
solutions ipn, let us check conditions on a parabolic boundary of the form [r,n] x {tq} U {f} x 
(''"Oi ''"i) U {n} X (tq, Ti), where ti > tq is a positive time to be chosen later. On the bottom we 
have that, thanks to the choice of C and (5.21 ), 

*(r, To) = C - ^ > 1 - Co > ij{r, tq) > Vn(r, tq) Vr € [f, n] , 

the last inequality following from the fact that V >Vn for r > 1. On the inner lateral boundary 
there holds 

^{f,T) = C-1-{t-to) = 1-^-1-{t-to) Vr G (To,ri) . 
r 2 r 

Now let us fix a small e > (let e < co/2). Assume that f is also larger than 2/e. Exploiting 
the local uniform convergence of ip to zero (Proposition 2.6), we know that if tq > ro(r,e) then 

e 



il^ir, t) < - Vr > Tq 



Hence, once r satisfies 



r < n = ro + 1 



Co 
2 



> ro + 1 - Co > ro 



it is guaranteed that 



^{f, r) > ^{f, n) = - > V'(r, r) > Vn(r, r) Vr G (rg, n) . 

On the outer lateral boundary {n} x (ro,ri) it is plain that ^'(n, r) > ij)n{n,T) = — oo. In fact 
this is the reason why we needed to suitably approximate V from below with the sequence Vn, 
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otherwise we would have not been able to compare ^ with ip in an outer lateral boundary. So, 
from the comparison principle, we get that ^' > -0^ in the region [f,n] x [ro,Ti]. In particular. 



i^ni^jTi) < ^'(r, Ti) < £ Vr G [f,n] . 
Passing to the limit in (5.22) as n — t- oo yields 

ipir, Ti) < ^(r, Ti) < £ Vr G [f, oo) 



Since ri = tq + 1 — co/2 — 1/f — e/2 and (5.23) holds for all tq > To(r,e), we have: 

'^{r, t) < e V(r, r) e [r, oo) x [f, oo) , 



(5.22) 
(5.23) 
(5.24) 



where r = ro(f,e) + l — co/2 — 1/f — e/2. Thanks to (5.24) and to the local uniform convergence 
of to zero, we conclude that 

limsup sup ip{r,T)<e. (5.25) 

T-!>oo re[0,oo) 



Because (5.25) is true for all e > 0, it holds for e = too and the proof is complete. □ 
The final step is to prove an analogous result for the relative error (j) = w'^/V — 1. Since 



the ideas are similar to the ones developed in the proofs of Lemma 5.3 and Proposition 5.4 
shall proceed in a concise way. 



we 



Lemma 5.5. Let tq > 0. Consider the function 

^r,T) = C-^-A{T-To) 
r 

If the positive parameters A, B, C, f satisfy the condition 
5^^e(iv-i)(i-e)(i,-i)f > (i + c')^-i ( - 

where e > is a small fixed number (e = 1/5 will work) and K^^ is a positive constant, then <1> 
is a supersolution to (5.2) in the region {(r, r) : <I>(r, r) > —1} Pi {(r, r) G (f j oo) x (ro,oo)}, 
provided r > f\/r^, where = r*(e, m, N) > 0. 



A 2(1 + C) 



m]T 



(5.26) 



(5.27) 



Proof. Condition (5.27) can be obtained reasoning as in the proof of Lemma 5.3 Indeed it is 
even easier to get to an inequality like (jslof with (1 - ^)V™-1 replaced by (1 + $)Vm-i and 
(1 — ^) replaced by (1 + $)^/™, because here we need not deal with the approximating sequence 
Vn- Then, in order to arrive at (5.27), one just uses the fact that by construction <I> < C □ 



Thanks to the (family of) barriers (5.26), we can now prove the analogue of Proposition 5.4 



for <j). 

Proposition 5.6. Let (p = jV — 1. There holds 



lim sup sup 1 

T-!>oo -re[0,oo) 



(r,r) < 0. 



Proof. Again, we proceed along the lines of the proof of Proposition 5.4 Fix A = B = 1, 
C = ci — 1/2, e < 1/2 and f large enough so that it satisfies f > f V r,,, V 2/e, where f complies 
with (5.27). These choices ensure that $ is a supersolution to (5.2 ) in the region (f , oo) x (tq, ti), 
with 

1 1 e 

Tl = Tn + Cl . 

U-r i 2 2 

Then, one takes tq = ro(f,e) so that (j){f,T) < e/2 for all r > ro(f,e); in this way $ and 
(j) are correctly ordered on [f,oo) x {tq} U {f} x {tq,ti). However, we have no control over 
their relation in an outer lateral boundary of the form {n} x (to,ti), for n > f. In order to 
overcome this difficulty it is enough to replace w by Wn < w, the latter being the sequence of 
rescaled solutions to (2.4), which vanish on {n} x (ro,ri). By construction (pn = w^!/ /V — 1 < (p 
and the equation solved by (pn is basically the same as (just replace (1 + (p)^/'^/T with 

(1 + (pn)^^'^/Tn), SO that for n large (that is, T„ close to T) $ is a supersolution also to such 
equation. Therefore from now on one proceeds exactly as in the proof of Proposition 5.4 □ 
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Going back to the original variables t and u{r,t), we have that Propositions 5.4 and 5.6 give 
convergence of u"^{r,t) to (1 — t/T)'"/(^~™')y(r) in relative error. From this fact the first main 
result stated in Theorem namely formula (1.5), clearly follows. 

6. Results for derivatives 

This section is devoted to proving the claimed results of Theorems |1.3| and |1.4| which 
deal with derivatives. In order to do that, we shall exploit a useful change of spatial variable 
and rescaling techniques in the spirit of [2^. 

Let 

' = '^'^ = 1 (sinho;)^-! 

so that 

as ar 

Notice that s ^ as r f oo, while s f oo as r 4- 0. Hence, the hyperbolic Laplacian of a regular 
function / = /(r) = /(r(s)) is 



A/ = (sinhr) 



-2(iV-l) 



ds2 • 



Let w be the solution to (2.3). Take tq > 1 and ro > R with R large enough (to be chosen 
later), and let sq = s(ro). Consider the rescaled function 

W{x,y) = (sinhro)^'w(r(so + x(sinhro)"^^"^)),To + y(sinhro)~(^~^H^-i)^ , (6.2) 

where (x,y) varies in the square (— a,a)^, a > being a fixed number to be set later. After 
straightforward computations one sees that W satisfies the following equation: 

(sinhro)2(^-i) 



(sinhr)2(^-i) 



(sinhro)-(^-iH^-i) 

+ — W . 

(l-m)T 



(6.3) 



a{r,ro) b{ro) 

From now on we shall fix a = a{N), independent of tq and tq, so small that s is forced not 
to go out of the interval {sq/2, 3sq/2) as x varies in the interval {—a, a). Such a choice of a is 
feasible. Indeed, notice first that from (6.1) one has 

lim 



C{N) > , 



(6.4) 



s{r) 
(sinhr)-(^-i) 

Moreover, s varies in the interval (sq — «(sinhro)^^^~^\ so + a(sinhro)~^'^~^-'). Hence, if ro > i? 
(namely, sq < s{R)) with R = R{N) large enough, a proper choice of a = a{N) (sufficiently 
small) ensures that s does not leave the interval (so/2,3so/2) as x S {—a, a). Still from (6.4) 
one deduces that 



lim r(s/2) - r{s) = C, lim r(3s/2) - r(s) = C 



(6.5) 



for suitable constants C and C. The choice of a above and (6.5) imply that, for tq = 'r(so) large 
enough (again, greater or equal than a given value R that depends only on A^), there holds 

1 



sinhro < sinhr < Ci sinhro Vr G (r(so/2) , r(3so/2)) 



(6.6) 



Ci being another positive constant. Gathering all this information, let us go back to equation 
(6.3). First, from the global Harnack principle (1.9) and from ( |6.6[ ), one deduces that 

\W{x,y)\<M V(x,y) G (-a,a)2 (6.7) 

for a positive constant M = M{uo,m, N) that does not depend on ro, tq. Moreover, notice that 

d (sinhr)^"^ d 

(sinhro)^"-*^ dr 



dx 



(6.8) 
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Hence, (6.3), (6.8) and the inequality 
1 



coshrg < coshr < Cicoshro Vr G (r(so/2), r(3so/2)) 



which can be proved as (|6.6|) above, give 
1 



Pk 



< 



d^a{r{x),rQ) 



<Pk, \b{ro)\<Q V(x,y) G (-a,a)^ V/c GNU {0}, 



(6.9) 



again for positive constants = Pfc(A^) and Q = Q{T{uo),m, N) that do not depend on rg, tq. 



The bounds ( |6.7| ) and (6.9) permit to conclude, as in the proofs of Theorem 1.1 and Lemmas 
3.1-3.2 of [21j (one uses the regularity results of |24| and |16| for bounded solutions to nonlinear 
parabolic equations like (6.3)), that the following estimates hold: 



(0,0) 



< M- 



l,k ' 



(0,0) 



< Ma fc VA; G N , 



(6.10) 



Mi^fc and M2^k being suitable positive constants depending only on M, a, T, m and A^. Go 
ing back to the original variables, we end up with the existence of positive constants Aij^ = 
Ai^kiuo, ni, N, e), A2^k = ^2,fc(^io, "m, N, e) such that 



d w 



Qj-k 



< A 



2,fce 



h,k e 



(6.11) 
(6.12) 



Vro > , Vro > e , VA; G N . 

In order to prove ( 6.1ip -(6.12 ) in the region tq > i?, tq > 1 one just uses (6.2), (6.8), (6.10) and 

the equality dr = (sinhro)~^^~^^('"~"^)dy recursively. To extend such bounds to tq G (0, i?) it 
is enough to apply the aforementioned regularity results to w directly, since in this region w 
is bounded away from zero. Finally, they also hold for tq > e up to admitting dependence of 
the constants vli^fc, ^2,fc on e as well through a and M (as a consequence of the global Harnack 
principle (1.9)). 

Using the fact that the solution V{r) to (1.4) and any of its derivatives V^^\r) behave like 
^-{N-i)r infinity (recall (1.3) and see formula (6.15) below), it is only a matter of tedious 
computations to check that estimates (6.11) imply that 



Qj.k 



<Aik Vr > e , VA; G N 



(6.13) 



(up to possibly different constants that we keep denoting as Ai^k), where ip is the relative error 
{p{r, t) = w{r, t) /y^/"^(r) — 1 (notice that it is different from the relative error (f) we dealt with 
in Section [5]). By exploiting the interpolation inequalities (we refer to |3H p. 130] or [9l App. 
3] for a short review on general interpolation inequalities) 



Mr)\\ 



i-k/j 



the bounds (6.13) and ( |1.5[ ) (that is, the fact that ||v5(t) 
also 



V/c, j G N : k < j , 

— )• as T — )• 00), we conclude that 



lim \<f{T)\ 



VA; G N . 



(6.14) 



Recalling (2.2), this is equivalent to formula (1.6) of Theorem [lT] As for estimates (1.10)- (1.11) 
of Theorem 1.3 just notice that they can be readily deduced from ( |6.11 )-(6.12 ) and again (2.2). 

Finally, let us investigate what (6.14) means in terms of spatial derivatives for w (or u). That 
is, we are going to prove formula (1.7) of Theorem |1.1[ First of all we show, by induction, that 

-a{N 



lim 



y"(r) 



1))^ Va / 0, VA; G N. 



(6.15) 
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Indeed, the existence of the hmit in (6.15) is a consequence of the chain rule and (1.3). Then, 



suppose that (6.15) holds for a given k. Consider the limit 



lim 



-(r) 



1 



dr 



lim 



1 



lim 



-(^) 



(^) ar^oo y°(r) y(r) -a(iV - 1) r^oo ^"(r) 



(6.16) 



Since the limit on the r.h.s. of (6.16) exists, by de I'Hopital's Theorem it coincides with 



lim 



this proves the inductive step. Because (6.15) is trivially valid for = 0, we conclude that it 
holds for all /c G N. 



Now recall that, for the k-th spatial derivative of c/?, one has the binomial formula 

d^^> fk^ 
J. 



Qj.k 



E 

j=0 



(6.17) 



where, to simplify notation, we use apexes to denote derivatives w.r.t. r. We are going to show 
that, for all /c G N, the ratio w^''\r,T)/V^^"^{r,T) converges in to a smooth function Gk{r) 
such that 

k 



lim Gfc(r) = (-1)* 



A^- 1 



m 



(6.18) 



This is clearly equivalent to proving (1.7). For k = (6.18) is exactly (1.5), with the choice 
Go = 1. For further derivatives we shall prove that (6.18) holds by induction. Indeed, suppose 
such result is true for all j < k — 1. From (6.14) and (6.17) we get 

^'^ (k-j) 



lim 

r^oo 



W 



ik)( 



Vr 



j=0 



(6.19) 



Thanks to the inductive step and (6.15), (6.19) implies 



lim 

r— >oo 



E 

j=0 



{k-j) 



So, we are left with proving that the function 

fc-i 



Gk{r) 



j=0 



G,(r) 



V- 



(k-j) 



V- 



m r 



complies with (6.18). From (6.15) and the inductive hypothesis this will follow provided there 
holds 

k 



N 



m 



an equality which is in fact true in view of the identity X]j=o 
Newton's binomial formula. 



1)'^ ^, valid by 



Remark 6.1. In order to obtain the bounds ( 1.14 )-( 1.15 ) of Theorem 1.4 (subcritical m), 
notice the quantities ci(t, •), C2(t, •) appearing in (1.13) can be taken to be continuous functions 
of t G (0,T), as a consequence of the method of proof of such inequality (see Sections |3] and |4]) . 
This fact is sufficient to prove, as above, that there exists a constant M (now depending on tq 
too) such that (6.7) holds true. This is enough in order to proceed along the previous lines and 
get (iri4l)-(irT5). 
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